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ABSTRACT 
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new  design  method  is  presented  which,  when  implemented,  uses  model  fol¬ 
lowing  and  full-state  feedback  to  keep  the  dominant  roots  of  a  system 
constant.  Under  favorable  circumstances,  it  can  do  this  even  in  the 
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presence  of  arbitrarily  large  parameter  uncertainties.  The  method  has 
the  attractive  feature  that  the  parameter-insensitivity  and  disturbance- 
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to  several  aircraft  flight  control  problems. 
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ABSTRACT 


The  problem  of  designing  multivariable  control  systems  is  addressed 
through  the  use  of  model -fol 1  owing  methods.  A  number  of  different 
model -fol 1  owing  techniques  are  discussed*  and  their  advantages  and  dis¬ 
advantages  are  presented.  It  is  shown  that  model -fol 1  owing  theories  can 
provide  useful  structural  insight  to  the  designer,  and  that  they  can 
help  to  integrate  the  methods  of  "classical"  and  "modern"  control.  A 
new  design  method  is  presented  which,  when  implemented,  uses  model  fol¬ 
lowing  and  full -state  feedback  to  keep  the  dominant  roots  of  a  system 
constant.  Under  favorable  circumstances,  it  can  do  this  even  in  the 
presence  of  arbitrarily  large  parameter  uncertainties.  The  method  has 
the  attractive  feature  that  the  parameter-insensitivity  and  disturbance- 
rejection  characteristics  of  the  system  can  be  selected  independently 
from  the  no-disturbance,  nominal -pi ant  performance.  Application  is  made 
to  several  aircraft  flight  control  problems. 
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Chapter  I 


INTRODUCTION 

1.1  BACKGROUND 

"Model  following"  refers  to  a  number  of  very  general  control -system 
design  techniques  which  can  be  used  to  solve  a  wide  variety  of  control- 
system  problems.  In  the  typical  model-following  problem,  it  is  desired 
to  force  some  or  all  of  the  states  of  a  plant  to  "match",  to  a  greater 
or  lesser  degree,  the  states  of  an  ideal  model.  This  model  could  be  an 
actual  hardware  system;  however,  it  is  normally  a  mathematical  model  of 
the  desired  closed-loop  plant  behavior. 

Two  of  the  most  common  applications  of  model  following  have  been  in 
the  area  of  aircraft  flight  control.  The  first  type  of  problem  involves 
the  design  of  the  flight  control  system  for  a  new  aircraft.  Frequently, 
by  referring  to  the  Military  Specifications  for  aircraft  flying  quali¬ 
ties,  through  simulator  studies,  or  from  past  experience,  the  designers 
have  established  desired  time-history  responses  to  pilot  commands.  The 
objective  is  then  to  select  the  control -system  structure  and  the  control 
gains  so  that  the  aircraft  will  have  responses  which  match  the  desired 
ones  as  closely  as  possible.  Such  a  problem  is  ideally  suited  to  the 
application  of  model-fol louing  techniques.  In  the  second  type  of  prob¬ 
lem,  we  would  like  to  make  one  kind  of  aircraft  "feel"  to  the  pilot  like 


a  different  kind  of  aircraft. 


There  are  several  reasons  why  we  might 


want  to  do  this.  For  example,  it  is  than  possible  to  conduct  systematic 
studies  of  different  aircraft  characteristics  to  determine  their  impact 


on  handling  qualities  or  the  usefulness  of  new  control  modes.  Also,  it 
is  often  very  helpful  to  be  able  to  provide  pilot  training  and/or  ini¬ 
tial  evaluations  for  a  neu  aircraft  before  the  aircraft  itself  is  avail¬ 
able.  By  using  special  "variable-stability  aircraft,"  preliminary  test¬ 
ing  has  been  accomplished  on  a  number  of  recent  development  programs, 
such  as  the  Supersonic  Transport  CSST),  the  F-16  fighter,  and  the  Space 
Shuttle,  long  before  the  final  prototypes  uere  completed.  Such  efforts 
are  classic  examples  of  model-fol lowing  problems  since  the  goal  is  to 
make  the  first  aircraft  (the  plant)  act  like  the  second  aircraft  (the 
model ) . 

Although  many  papers  have  been  written  on  the  model -fol I  owing  prob¬ 
lem,  the  question  of  how  model  following  relates  to  control-system  sen¬ 
sitivity  has  not  been  fully  addressed.  The  purpose  of  this  thesis  is  to 
review  the  various  types  of  model-following  techniques,  to  show  how  they 
are  related  to  each  other  and  to  other  design  methods,  and  to  show  how 
model  following  can  be  used  for  designing  insensitive  controllers. 

1.2  LITERATURE  REV  I  Ell 

There  are  two  ways  to  approach  the  actual  implementation  of  a  model- 
fol  louing  system.  In  the  first  method,  Implicit  Model  Following  (also 
known  as  "Model  in  the  Performance  Index"),  the  model  equations  are  only 
used  during  the  design  process  to  aid  in  the  determination  of  the  con¬ 
trol  gains.  After  the  design  is  complete,  the  model  is  no  lonaor 


needed,  and  it  is  not  necessary  to  implement  the  model  equations  in  the 
actual  system.  The  second  method,  knoun  as  Explicit  Model  Following 
(also  called  "Real  Model  Following"  or  "Model  in  the  System"),  requires 
the  model  equations  to  be  incorporated  as  part  of  the  control  system. 
The  model  states  are  then  assumed  to  be  available,  along  with  some  or 
all  of  the  plant  states,  for  formulation  of  the  control  signals. 
Kreindlerl 161  and  others  have  compared  the  explicit  and  implicit  model- 
following  methods  for  particular  aircraft  examples;  however,  conclusions 
on  the  overall  merits  of  the  two  approaches  have  usually  been  rather 
subjective. 

Uhether  a  system  uses  an  implicit  model -fol 1  owing  approach  or  an  ex¬ 
plicit  one,  it  is  clear  that  the  ideal  result  would  be  to  have  the  plant 
dynamics  exactly  match  the  model  dynamics.  Such  a  capability  is  re¬ 
ferred  to  as  "perfect  model  following."  If  perfect  model  following  is 
possible  (and  the  control  gains  are  selected  correctly),  then  uhen  the 
plant  and  model  are  started  at  the  same  initial  conditions  and  with  no 
unknown  disturbances,  one  or  more  of  the  plant  states  will  identically 
match  the  corresponding  model  states.  Erzbergerl 10 1  developed  suffi¬ 
cient  conditions  and  an  expression  for  the  control  gains  to  obtain  per¬ 
fect  model  following  for  one  class  of  implicit  model-following  systems, 
while  Chanl5]  and  Morrow  and  Balasubramanian[20]  gave  the  appropriate 
control  laws  for  the  explicit  formulation.  Asseo[31  and  Motyka  and 
Rynaski[221  extended  some  of  these  perfect-model-fol lowing  ideas,  ulnle 
CurranlSl  developed  the  related  concept  of  "equicontrol labil ity",  and 
showed  how  perfect-model-following  theories  could  be  made  more  useful  by 
allowing  slight  modifications  of  the  model. 


Ty!er[291  examined  the  use  of  optimal  control  theory  for  solving 
model-fo! lowing  problems.  He  was  followed  by  a  number  of  researchers, 
including  Kriendlerl  161,  Gran,  Berman,  and  Rossillll,  Kriechbaum  and 
StinemanI 17],  Marklandl 19],  Tiroshi  and  E!1iot[27],  and  Yore[32]. 
Rediess  and  Uhitaker[25]  and  Peled[24l  also  used  models  to  formulate 
their  optimal -control  cost  functions,  but  their  methods  were  really  par¬ 
ameter-optimization  schemes  which  required  the  designer  to  determine  the 
control  system's  structure  and  the  appropriate  form  of  any  required  com¬ 
pensation.  Trankle  and  Brysonf28]  developed  a  variation  of  earlier 
methods  which  they  found  useful  for  solving  some  classes  of  model-fol¬ 
lowing  problems. 

Many  of  the  above  authors  stated  their  "intuitive"  feelings  that  ex¬ 
plicit  model-following  systems  had  potential  for  reducing  the  sensitiv¬ 
ity  of  control  systems  to  parameter  uncertainty.  Nevertheless,  there 
have  been  feu  if  any  attempts  to  quantify  this  capability,  or  to  formu¬ 
late  an  appropriate  design  method.  Uinsor  and  Roy] 31 1  attempted  to 
tackle  the  sensitivity  issue,  but  their  design  method  was  complex  and 
was  limited  to  small  parameter  variations.  In  addition,  they  concluded 
that  it  was  not  possible  to  desensitize  more  than  one  variable  at  a 
time.  Landau  and  CourtiolllS]  apparently  endorsed  Winsor  and  Roy's 
findings  and  stated  their  belief  that  an  adaptive  system,  even  with  its 
inherent  complexity  and  nonlinearities,  was  the  best  approach. 
Kamiya[15]  was  probably  the  first  to  study  the  use  of  mode  1 -f o  1 1  oumg 
systems  specifically  for  parameter-insensitive  control;  however,  he  used 
a  different  structure  than  the  one  proposed  in  this  work. 
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Since  one  of  the  major  contributions  of  this  uork  relates  to  parame¬ 
ter-insensitive  controllers,  it  is  important  to  mention  briefly  some  of 
the  other,  non-model -fol lowing  efforts  in  this  area.  Palsson  and 
Uhitaker[23]  and  Hadass[121  developed  the  technique  of  treating  the  un¬ 
certain  parameters  as  random  vectors.  Harvey  and  Pope! 13]  and 
Vink!er[30]  compared  several  different  methods  in  their  works,  while 
Shenkar[26]  and  Ashkenazi [21  extended  the  most  promising  methods  for  the 
output-feedback  case.  In  general,  all  of  the  methods  proposed  required 
large,  time-consuming  computer  programs  and  used  gradient-search  algor¬ 
ithms  to  obtain  the  final  results.  Even  so,  the  resulting  systems  had 
the  unfortunate  characteristic  that  the  system  performance  and/or  con¬ 
trol  use  for  the  nominal  system  was  typically  degraded  from  what  would 
have  been  obtained  had  the  design  been  accomplished  without  regard  to 
the  parameter  variations. 

1.3  CHAPTER  OUTLINES 

Chapter  II  defines  perfect  model  following  and  discusses  how  it  re¬ 
lates  to  other  methods  such  as  pole  placement.  The  criteria  for  achiev¬ 
ing  perfect  model  following  are  presented,  and  it  is  shown  that  they  can 
provide  useful  insight  into  the  structure  of  a  control -system  problem, 
even  when  other  design  methods  are  used. 

Chapter  III  shows  how  we  can  combine  model-following  techniques  and 
optimal  control  theory  to  take  advantage  of  the  respective  strengths  of 
"classical"  and  "modern"  control.  Both  state-matching  and  dynamics¬ 
matching  systems  are  developed,  and  the  advantages  and  disadvantages  of 
each  type  of  system  are  presented. 
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Chapter  IV  begin*  with  *  review  of  peat  efforts  in  parameter-inaensi 
tive  control.  A  new  design  method  using  model  following  is  then  pre 
sented  which  is  shown  to  have  significant  advantages  over  previous  tech 
niques. 

Chapter  V  gives  conclusions  and  recommendations  for  future  research. 

1.4  SUMMARY  flF  CONTRIBUTIONS 

1.  A  number  of  different  model-following  design  techniques  are 
evaluated,  and  their  advantages  and  disadvantages  and  rela¬ 
tionships  to  other  methods  are  presented. 

2.  The  relationship  between  feedback  and  feedforward  control  is 
demonstrated. 

3.  A  neu  design  method  is  presented  which  uses  model  following 
and  full-state  feedback  to  keep  the  dominant  roots  of  a  system 
constant.  Under  favorable  circumstances,  it  can  do  this  even 
in  the  presence  of  arbitrarily  large  parameter  uncertainties. 
The  method  has  the  attractive  feature  that  the  parameter-in¬ 
sensitivity  and  disturbance-rejection  characteristics  of  the 
system  may  be  selected  independently  from  the  no-disturbance, 
nominal -parameter  performance.  The  method  is  applicable  to 
systems  in  which  the  plant  can  be  made  "well -damped"  and  "sub¬ 
stantially  faster"  than  the  final  system  characteristics  re¬ 
quired  by  the  model . 
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Chapter  II 

PERFECT  MODEL  FOLLOWING  AS  A  MULTIVARIABLE  DESIGN  TOOL 

2.1  REVIEW  OF  POLE-PLACEMENT  TECHNIQUE 

The  use  of  full-state  feedback,  either  by  measurement  or  by  estima¬ 
tion,  gives  the  control -system  designer  great  flexibility,  since  the 
roots  of  the  closed-loop  system's  characteristic  equation  may  be  placed 
arbitrarily  with  the  proper  choice  of  control  gains.  This  has  led  to  a 
very  simple,  easy-to-use  design  technique  for  single-input,  single-out- 
put  systems.  In  this  technique,  known  as  pole  placement,  the  character¬ 
istic  equation  of  the  desired  closed-loop  system  is  compared  to  the 
characteristic  equation  of  the  open-loop  plant,  modified  so  that  there 
are  feedback  gains  on  all  of  the  states.  By  matching  the  coefficients 
of  the  Laplace  variable  "s",  a  set  of  algebraic  equations  results  which 

i 

allows  for  solution  of  the  control  gains. 

Let  us  take  a  simple,  second-order  system  as  an  example.  The  differ¬ 
ential  equation  of  motion  of  an  open-loop  plant  is  given  by  x+x=u,  or  in 


state-space  form, 


Suppose  the  desired  system  has  a  natural  frequency  of  2  rad/sec  and  a 
damping  ratio  of  .7.  The  desired  characteristic  equation  is  then 

s2  +  2 . 8s  +  4  =  0 . 


Using  full-state  feedback,  the  control  lau  is  given  by  u=Cx,  so  ue  have 


■ 

' 

* 

■ 

xi 

0 

1 

Xl 

X2 

-  1+Ci 

C2 

. 

X  2 

with  a  characteristic  equation  given  by 


s2  -  c 2 s  +  1  -  c i  =  0 . 


Comparing  coefficients,  ue  see  that 

— c  2  =  2.8 

and 


1  -  Cl 


Cl  =  -3 


and 


c2  =  -2.8. 


Unfortunately,  the  technique  breaks  down  for  multiple  inputs  and  out¬ 
puts.  In  the  above  example,  there  are  tuo  degrees  of  freedom  from  the  ! 

two  control  gains,  and  tuo  equations  resulting  from  comparing  the  coef¬ 
ficients  of  "s".  Houever,  with  additional  controls,  additional  degrees 
of  freedom  are  added  through  additional  control  gains,  but  there  are  no 
additional  equations,  thus  preventing  a  unique  solution.  Suggestions 
for  circumventing  this  problem  include  arbitrarily  setting  some  of  the 
control  gains  equal  to  2ero,  or  assuming  that  some  of  the  entries  in  the 
dynamics  matrix  are  small  enough  to  be  neglected.  In  either  case,  sev¬ 
eral  iterations  are  usually  required,  since  there  is  no  guarantee  that 
the  resulting  system  will  have  desirable  response  characteristics. 

i 

>■ 

2.2  DERIVATION  flF  PERFECT- MO DEL-FOLLOWING  CONTROL 

'  Suppose  that  rather  than  matching  the  characteristic  equations,  ue 

I  try  to  match  the  entire  dynamics  matrices.  Repeating  Ersbergor's 

derivation! 10 ] ,  but  with  our  notation,  ue  let  the  plant  be  described  by 

|  x  =  Fx  +  Gu 

[  y  =  Hx. 


The  model,  which  describes  the  desired  system  behavior,  is  given  by 


=  F  m* 


for  y  and  xa  to  have  the  same  dynamics,  ue  require  that 


y  =  f  my  =  r„,Hx. 


Also, 


y  =  Hx  =  HFx  +  HCu. 


Thus,  ue  have 


HFx  +  HGu  =  FmHx. 


Solving  for  u,  we  obtain 


u  =  I(HG)t(HG)r,(HG)t(raH>HF)x. 


Ue  define  the  pseudo-inverse  of  HG  as: 


-  10  - 


(HG  )*  =  [(HG)t(HG)l',(HG)t. 


Then , 


U  =  (HG)*(FmH  -  HF)x. 

By  substituting  this  expression  for  u  into  the  equation  for  y,  ue  get 

y  =  HFx  +  (HG) (HG)*(FmH-HF)x. 

Then,  in  order  for  y  to  equal  Fmy,  ue  must  have 

HFx  +  (HG)(HG)*(FmH  -  HF)x  =  FmHx 
l  HF  +  (HG)  (HG)*'.F„H-HF)  -  FmH]x  =  0 
(HG)(HG)*(FmH-HF)  -  I ( F^H  -  HF)  =  0 
[ ( HG ) ( HG ) *  -  I ] [ F mH  -  HF]  =  0. 


The  final  equation  above,  known  as  Erzberger's  condition,  gives  a  suffi¬ 
cient  condition  for  exactly  matching  the  dynamics  of  the  plant  and  the 
model  for  this  class  of  problem.  Such  matching  is  known  as  "perfect 
model  following". 
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14 e  can  generalize  the  above  procedura  to  handle  inputs  to  the  model, 
as  well  as  the  matching  of  only  a  subset  of  model  states.  Let  the  plant 
again  be  described  as 


x  =  Fx  +  Gu 

y  =  Hx . 


The  desired  system  dynamics  are  given  by  the  model  equations: 


X  m  -  F  »X  m  +  G  nU  m 

y.  =  HmX*. 


We  would  like  to  have 


y  =  y*- 


Thus, 


HFX  +  HGU  =  IIbF»Xi«  +  HmGmUm 
HGu  •  H mF  fflX in  *  HFx  ^  HmGmUm 
U  =  ( H G ) * H mF mx m  -  CHGJ-HFX  +  (HG)*H„Gmum. 
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Or 


U  =  C!X  +  C2X„  +  C  3U  n  p 

where 

Ci  =  -(HG)*HF 
C2  =  CHG)*HnF« 

and 

C3  =  (HG)*Ha,Gn. 

Notice  that  this  general  expression  for  the  perfect-model -fol 1  owing  con¬ 
trol  contains  three  terms:  a  feedback  of  plant  states  x,  a  feedforward 

of  model  states  xm,  and  a  feedforward  of  model  inputs  um.  However,  this 
does  not  mean  that  we  are  restricted  to  this  type  of  structure.  If  we 
are  able  to  successfully  make  a  plant  state  exactly  follow  a  model 
state,  then  obviously  the  two  signals  are  equivalent,  and  we  can  substi¬ 
tute  one  of  them  for  the  other  in  the  control  law.  For  example,  suppose 
we  are  working  a  problem  in  which  the  number  of  states  to  be  matched  is 
equal  to  the  number  of  model  states  (i.e.,  Hm  is  square).  Then  Hr,'1 
wi l I  normal ly  exist. 

Since  Hx  =  Hmxm, 


The  perfect-model-fol lowing  control  lau  can  then  be  written  as 

U  =  <HG)*Hn,F »Hm‘ 1  Hx  -  (HG)*HFx  +  (HG)*H.G«u« 

=  (HG )*[ HmFBHn*  *  H  -  HFjx  +  (HG)*H.G,u„. 

Thus  for  this  case,  ue  can  achieve  perfect  model  following  with  only 
feedback  from  the  plant  states,  without  having  to  implement  the  model 
system.  For  Hm=I  and  Gm=0.  ue  are  left  with 

u  =  (HG)*[F*H  -  HFlx, 

which  is  the  control  lau  derived  ,by  Erzberger. 

If  the  number  of  states  to  be  matched  is  equal  to  the  number  of  plant 
states  (i.e.,  H  is  square),  than  H’ *  will  normally  exist. 

Thus, 

x  =  H-'Hwx„ 

and 


U  =  (HG)*ilnFmx„  -  CHG)*HFH*1Hmxra  +  (H0)*HmGmum 


=  (HG)*[HaFa  -  HFH* ’  Ha]xa  ♦  (HG)*HaGaua. 


This  control  lau  can  therefore  be  implemented  with  only  feedforward  from 
the  model  states,  without  the  need  for  measuring  any  plant  states. 

Finally,  if  al 1  of  the  plant  states  are  matched  to  al 1  of  the  model 
states,  then  both  H  and  H„  are  square  and  invertible,  and  either  of  the 
above  control  laws  can  be  used.  As  a  result,  it  is  possible  to  imple¬ 
ment  the  control  system  with  either  feedforward  or  feedback,  or  a  combi¬ 
nation  of  the  two. 

By  substituting  the  general  expression  for  u  into  the  equation  for  y, 
we  get 


y  =  HFx  +  (HG) (HG )* [  H«F  nxm  -  HFxl  +  (HG)  (HG)*HmGaUm. 


We  would  also  like  to  have 


e 

y  =  HraFmXm  +  HmGnUin. 


Equating  the  two  expressions. 


HFx  +  (HG ) (HG)*l HmFmXm  -  HFx]  +  (HG) (HG)*HmGmUa  =  HmFaxm  +  HmGmu 

(HG) (HG)*[HmF  mxm  -  HFxl  +  [HFx  -  HmF„xm] 

+  (HG)(HG)*HaGmUm  -  HmGmU,,,  =  0 
-  15  - 


( (HO  (HO* 


1  1 1  HaFaXa  -  HFx  )  +  UHGHHGJ-  - 


I  ]HaGau 


=  0. 


Since  ua  is  arbitrary,  ue  must  require  that 

t(HGHHC)*  -  1 1 1 HaFaxa  -  HFxJ  =  0 

and 

ICHGXHG)*  -  I]HaCa  =  0. 

The  above  conditions  give  the  general  criteria  for  perfect  model  follow¬ 
ing  to  exist.  By  examining  specific  cases,  it  is  possible  to  use  these 
relations  to  get  meaningful  conditions.  For  example,  if  H„=l  and  Ga=0, 
ue  get 

[(HGHHO*  -  !l[Faxa  -  HFx  1  =  0. 

Since  xa=Hn* 'Hx=Hx,  ue  can  write 

[  (HG)  (HO*  -  I  1 1  F mH  -  HF  lx  =  0. 

Therefore, 

l  (HG)  (HO*  -  I  ]  I  F  mH  -  HF]  =  0, 
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which  is  Erzberger's  condition. 


2.3  INTERPRETATION  OF  CONDITIONS  FOR  PERFECT  MODEL  FOLLOWING 

It  is  interesting  to  note  that  a  sufficient  condition  for  perfect 
model  following  to  exist  is 

l(HGHHG)*  -11=0. 


First,  let  us  define  the  term  "directly  controllable"  to  mean  that  a 
state  has  at  least  one  non-zero  entry  in  the  corresponding  row  of  the 
control  matrix,  G.  For  example,  in  the  system 


Xl 

xz 


Ft  i 
Fzi 


x 2  is  directly  controllable,  while  xi  is  not.  The  fact  that  a  state  is 
directly  controllable  therefore  means  that  it  is  possible  to  modify  that 
row  of  the  dynamics  matrix.  With  this  definition,  the  above  condition 
can  be  understood  as  requiring  that  each  plant  state  to  be  matched  must 
be  directly  controllable,  and  that  the  number  of  controls  be  equal  to 
the  number  of  plant  states  to  be  matched.  The  reason  that  this  is  a 
sufficient  condition  and  not  a  necessary  one  relates  to  the  fact  that 
even  if  a  state  is  not  directly  controllable,  sometimes  the  row  of  the 


plant  dynamics  matrix  is  already  the  same  as  the  corresponding  row  in 
the  model  dynamics  matrix,  so  no  additional  control  effort  is  required 
to  force  the  match. 

Even  though  the  perfect-model -fol 1  owing  gains  are  straightf oruard  and 
easy  to  calculate,  it  might  at  first  appear  that  the  criteria  to  be  sa¬ 
tisfied  are  quite  stringent,  and  that  perfect  model  following  is  there¬ 
fore  rather  limited  in  use.  However,  with  some  slight  modifications, 
the  theory  can  be  extended  to  cover  a  much  broader  range  of  problems. 

To  restate  the  sufficient  condition  for  perfect  model  following: 


Plant  states  can  be  made  to  match  model  states  exactly 
if  they  are  directly  controllable  and  if  the  number  .of 
controls  is  equal  to  the  number  of  states  to  be  matched. 


If  the  number  of  controls  is  greater  than  the  number  of  states  to  be 
matched,  one  can  simply  delete  the  extra  controls.  If  the  number  of 
controls  is  less  than  the  number  of  states  to  be  matched,  then  the 
mathematics  does  yield  a  solution,  but  it  represents  only  an  approximate 
match  of  the  plant  and  model  states.  Actually,  it  is  a  "least-squares" 
compromise  of  the  desired  matches.  Whether  this  represents  an  accepta¬ 
ble  solution  depends  upon  how  different  the  open-loop  plant  and  model 
dynamics  are. 
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In  many  aircraft,  actuators  ara  used  to  drive  the  control  surfaces. 


For  such  systems,  the  primary  plant  states  are  not  "directly  controlla¬ 
ble"  as  defined  above,  thus  precluding  the  straightforward  application 
of  perfect-model-fol lowing  theory.  Fortunately  however,  the  dynamics  of 
aircraft  control  actuators  are  usually  much  faster  than  the  vehicle  dy¬ 
namics.  As  a  result,  they  can  normally  be  neglected  for  the  initial 
control -system  designs. 

Even  if  a  plant  state  is  not  directly  controllable,  it  will  still 
match  a  model  state  if  its  dynamics  are  the  same  as  the  model  state,  and 
the  state(s)  through  which  it  is  indirectly  controlled  are  made  to  match 
the  model.  For  example,  in  aircraft  equations  of  motion,  roll  angle  is 
normally  just  the  integral  of  roll  rate  for  both  the  plant  and  the 
model.  As  a  result,  even  though  roll  angle  is  not  a  directly  controlla¬ 
ble  state,  it  will  match  the  model  if  roll  rate  is  made  to  match  (assum¬ 
ing  the  initial  conditions  are  the  same).  If  we  have  a  situation  where 
a  plant  state  is  not  directly  controllable  and  there  are  no  correspond¬ 
ing  higher-order  states  in  the  model,  then  it  will  be  necessary  to  can¬ 
cel  the  extra  plant  states  in  order  to  get  an  exact  match.  For  example, 
the  only  way  to  get  a  1/s2  system  to  behave  like  a  1/s  model  is  to  some¬ 
how  get  rid  of  an  "s",  such  as  by  adding  a  differentiator.  With  a  noisy 
system,  a  solution  of  this  type  would  probably  not  be  practical.  Alter¬ 
natively,  as  Curran  has  shown[8l,  the  plant  can  be  made  to  match  a  modi¬ 
fied  model,  which  consists  of  the  original  model  plus  a  number  of  addi¬ 
tional  arbitrarily  fast  roots.  In  other  words,  we  could  have  our  1/s2 
system  behave  like  a  model  with  dynamics  given  by  1OOO/s(s+1O0O) .  Of 
course,  the  faster  the  additional  roots  are  (and  hence  the  closer  the 
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modified  model  is  to  the  original  one)*  the  greater  the  control  gains 
and  control  effort  which  will  be  required. 

Even  if  the  actual  selection  of  control  gains  in  a  system  is  done  by 
a  different  method*  applying  perfect-model-fol lowing  philosophy  can  of¬ 
ten  allow  the  control -system  designer  to  gain  some  helpful  structural 
insight  into  the  problem.  For  example,  the  number  of  controls  should  be 
based  on  the  number  of  states  for  which  a  particular  time  response  is 
required  (as  opposed  to  simply  being  stable).  Even  though  the  standard 
control -theory  criteria  may  indicate  that  all  of  the  states  of  a  system 
are  "controllable"  with  only  a  single  control,  if  it  is  desired  to  spe¬ 
cify  precisely  the  responses  of  more  than  one  state,  additional  controls 
will  probably  be  necessary.  Using  an  aircraft  as  an  example,  if  the 
elevator  is  the  only  longitudinal  control,  then  it  is  only  possible  to 
guarantee  a  particular  response  for  a  single  output  variable  -  whether 
it  be  pitch  rate,  or  normal  acceleration,  or  a  linear  combination  of  the 
two.  If  it  is  desired  to  specify  the  responses  for  both  pitch  rate  and 
normal  acceleration,  then  a  second  longitudinal  control,  such  as  a  ca¬ 
nard  or  maneuver  flaps,  would  be  required.  Along  the  same  lines,  the 
contro I -system  designer  should  strive  to  have  each  state  of  interest  be 
directly  controllable.  If  this  is  not  possible  because  an  actuator  is 
in  the  system,  the  control-system  designer  should  push  for  as  "fast"  an 
actuator  as  possible.  This  not  only  makes  the  designer's  job  easier,  it 
also  allows  him  to  design  a  control  system  which  will  provide  the  maxi¬ 
mum  flexibility  and  best  total  performance.  In  the  end  of  course,  ccst 
or  other  constraints  may  not  permit  the  control -system  designer  to  have 
the  final  say. 
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2.4  EXAMPLE :  LATERAL-DIRECTIONAL  DECOUPLING  OF  A  T-33  AIRCRAFT 

As  an  example  of  the  application  of  perfect-model-fol lowing  theory, 
consider  the  problem  of  decoupling  the  1 ateral -di rectional  axes  of  a 
T-33  aircraft.  This  example  has  been  solved  by  geometric  methods  in  [6] 
and  with  an  explicit  model -fol 1  owing  technique  in  [20]. 

Lateral-directional  decoupling  refers  to  the  desire  to  make  the  yau 
and  sideslip  motions  of  an  aircraft  independent  from  its  rolling  mo¬ 
tions.  If  this  can  be  accomplished,  the  pilot  workload  is  significantly 
reduced,  resulting  in  improved  safety  for  commercial  operations  or 
greater  ueapons-del i very  accuracy  for  military  aircraft.  The  particular 
aircraft  for  this  example  is  a  T-33  which  has  been  modified  for  flying- 
qualities  research.  In  addition  to  the  normal  rudder  and  aileron  con¬ 
trol  surfaces,  it  has  hydraulically  operated  "drag  petals"  on  the  wing- 
tip  fuel  tanks.  These  surfaces  can  be  differentially  extended  in  flight 
in  order  to  generate  significant  yawing  moments,  if  required. 

The  equations  of  motion  of  the  unaugmented  T-33  are  given  by 


x  -  Fx  *  Gu, 


where 


x*  =  (p  *  r  3] 


-3.18  0 


.63 


-10.6 


1. 

0. 

0. 

0. 

-.06 

0. 

-.27 

4.  18 

.022 

.0644 

-.998 

-.151 

-14.4 

1.5 

0. 

0. 

0. 

0. 

0. 

-2.59 

-.96 

0. 

.037 

0. 

and 


ut  =  ISa  6r  Spl. 


The  state  vector  consists  of  roll  rate,  roll  angle,  yaw  rate,  and  side¬ 
slip  angle,  while  the  control  vector  is  made  up  of  aileron  position, 
rudder  position,  and  drag  petal  position. 

The  model  equations  were  obtained  by  replacing  the  undesirable  cou¬ 
pling  terms  in  the  plant  equations  with  zeros.  The  resulting  model  is 
described  by 


Xffl  -  f  »*m  +  UnUju 


\ 

f 


r 


"I 


Fm  = 


Cm  = 

and 

u„* 

Since  there  are  three  controls,  and  the  roll  rate.  yaw  rate,  and 
sideslip  angle  are  all  "directly  controllable,"  we  are  guaranteed  that 
all  three  of  these  states  can  be  made  to  match  the  model  states  exactly. 
In  addition,  since  roll  angle  is  simply  the  integral  of  roll  rate,  it 
will  also  match  if  roll  rate  is  made  to  match.  Thus,  we  are  able  to 
match  exactly  all  four  plant  states  to  the  four  model  states.  Note  that 
this  system  may  be  implemented  either  with  feedforward  or  with  feedback. 
The  plant-state-feedback  or  model -state-feedforward  gains  are 
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The  gains  on  the  model  inputs  are 


1. 

.  1042 

0. 

0. 

1. 

0. 

0. 

-2.6979 

1. 

As  expected,  the  closed-loop  dynamics  matrix  using  these  gains  ex¬ 
actly  matches  the  model  dynamics  matrix. 


2.5  EXAMPLE: 


CONTROL  OF  A  B-26  AIRCRAFT 


In  the  previous  example,  the  plant  had  three  controls  and  three  di¬ 
rect  1  y-control  1  abl  e  states,  uhich  made  it  possible  to  match  exactly  all 
of  the  model  states.  Houever,  often  the  designer  will  be  faced  with  a 
problem  uith  feuer  controls  than  states.  Such  a  case  is  examined  here. 
The  plant  equations  represent  the  lateral-directional  equations  of  mo¬ 
tion  of  a  B-26  aircraft: 

x  =  Fx  +  Gu, 


x  *  =  [?  p  B  r ] 
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The  state  vector  consists  of  roll  angle,  roll  rate,  sideslip  angle,  and 
yaw  rate-,  while  the  control  vector  is  made  up  of  rudder  position  and 
aileron  position. 

The  model  was  determined  from  simulator  studies  to  provide  improved 
aircraft  handling  qualities: 


f  m  ^  G  m  * 


ui  th 


=  Pm  Bm  r  m  1 
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1. 

0. 

0. 

0. 

-1. 

-73.  14 

3.  18 

.086 

0. 

-.11 

-1. 

.0086 

.086 

8.95 

-.49 

G«  =  G 

and 

u*  *  =  (6rn  Sami . 


If  we  try  to  apply  Erzberger's  condition  for  perfect  model  following, 
ue  find  that  it  is  not  satisfied;  thus  both  Erzberger ( 10 ]  and  Chan(51 
concluded  that  "perfect  model  following  is  not  possible."  However,  both 
authors  proceeded  to  solve  for  the  perfect-model -fol lowing  gains,  and 
found  that  since  the  plant  and  model  equations  were  not  too  different,  a 
reasonably  close  "compromise"  match  of  the  four  states  could  be 
achieved . 

The  gains  for  this  case  were  determined  to  be 


Ci  or  C2 


-.0034 

0. 


-.1110 

-.4930 


-.3706 

17.4911 


-.0846 

-1.0128 


Since  the  plant  and  model  control  matrices  are  the  same,  the  gain  ma¬ 
trix  from  the  model  inputs,  C3,  is  simply  the  identity  matrix. 
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In  this  particular  problem,  the  plant  and  model  equations  uere  quite 
similar,  and  the  compromise  solution  which  attempted  to  match  all  four 
states  is  probably  acceptable.  Obviously,  this  will  not  always  be  the 
case.  Therefore,  it  is  important  to  note  that  by  using  the  generalized 
perfect-model-fol lowing  control  law,  it  Is  possible  to  achieve  perfect 
model  following  for  a  subset  of  the  plant  states.  Since  there  are  two 
controls,  we  are  guaranteed  to  be  able  to  match  two  states;  houever, 
since  roll  angle  is  simply  the  integral  of  roll  rate,  it  will  be  matched 
if  roll  rate  is  matched.  If  we  want  to  match  roll  angle,  roll  rate,  and 
yaw  rate,  the  selector  matrix  becomes 

10  0  0 
H  =  0  1  0  0 

0  0  0  1 

The  resulting  feedback  gain  matrix  for  the  plant  states  is 

0.  -.1084  .8749  -.1786 

Ci  = 

0.  -.7494  -1.2148  -.1995 

B 

while  the  feedforward  gain  matrix  for  the  model  states  is 


c  2  = 


0034 


-.0027  -1.2455  .0940 


0.  .2558  18.7059  -.8133 


Roll-angle  time  histories  of  the  4-state  "compromise"  match  and  the 
3-state  "exact"  match  are  given  in  Figure  1.  The  results  for  roll  rate 
and  yau  rate  are  similar,  with  the  3-state-match  design  duplicating  the 
model  system  exactly.  Although  ue  do  not  have  any  guarantees  as  to  uhat 
kind  of  response  ue  uill  get  for  the  state  ue  did  not  try  to  match 
(sideslip  angle),  in  this  particular  problem  it  is  very  close  to  the 
model  response. 

Obviousfy,  other  choices  of  H  are  possible  if  it  is  desired  to  match 
a  different  subset  of  the  states.  Essentially,  this  control  lau  uses 
the  feedback  from  the  plant  states  to  cancel  the  plant  dynamics  and  the 
feedforuard  from  the  model  states  to  obtain  the  desired  responses.  Note 
that  using  Erzberger's  control  lau  uith  this  type  of  problem  uill  not 
yield  exact  matches,  even  uith  the  same  selector  matrix  as  the  one  used 
here . 
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Figure  Is  Comparison  of  B-26  Roll-Angle  Responses 
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Chapter  III 


NO  DEL  FOLLOWING  AND  OPTIMAL  CONTROL 

3.  1  SELECTION  fl£  WEIGHTING  nATRICES  £fl&  OPT  1  UAL  £gNJML  P.RQB1EC!S 

Many  control -system  problems  can  be  formulated  in  such  a  uay  that  it 
is  desired  to  find  the  control  gains  which  will  minimize  a  cost  function 
of  the  form 


<9 


J 


J (x*Ax  +  u*Bu)dt. 
0 


where  A  and  B  are  weighting  matrices  which  penalize  state  excursions  and 
control  usage  respectively.  The  mathematics  for  the  solution  of  such 
"optimal  control"  problems  is  well  known,  and  is  easily  accomplished  on 
a  digital  computer.  The  resulting  control  laws  have  a  number  of  favora¬ 
ble  properties,  such  as  guaranteed  stability.  In  addition,  because  of 
the  matrix  structure,  multivariable  problems  are  no  more  difficult  to 
solve  than  single-input,  single-output  ones.  Nevertheless,  optimal  con¬ 
trol  techniques  have  not  seen  widespread  acceptance  in  industry,  parti¬ 
cularly  for  aircraft  control  problems.  Part  of  the  reason  for  this  may 
lie  in  the  relative  neuness  of  the  optimal  control  methods,  and  in  the 
availability  of  skilled  classical  designers  whose  experience  and  insight 
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allow  them  to  obtain  excellent  results  with  Root-Locus  or  Bode-Plot 
techniques,  even  for  very  complicated  problems.  However,  some  of  the 
reluctance  to  use  the  modern  control  methods  may  be  due  to  the  diffi¬ 
culty  of  selecting  the  weighting  matrices  A  and  B  in  order  to  obtain  the 
desired  results.  Increasing  the  weighting  on  a  state  tends  to  speed  up 
the  response  of  that  state,  while  increasing  the  weighting  on  a  control 
tends  to  cut  down  on  the  use  of  that  control.  However,  it  is  sometimes 
difficult  to  translate  such  specifications  as  frequency,  damping,  rise 
time,  or  overshoot  into  appropriate  choices  for  the  weighting  matrices. 
Thus,  many  design  iterations  may  be  required  before  acceptable  results 
are  obtained,  if  the  system  requirements  are  ones  of  this  type.  As  with 
classical  design  techniques,  a  great  deal  of  experience  and  insight  is 
found  to  be  helpful,  if  not  required,  to  achieve  the  desired  results. 

In  an  attempt  to  provide  some  guidance  for  the  selectTBn  of  the 
weighting  matrices,  Bryson's  Rule  was  formul atedl 4] .  Bryson's  Rule 
states  that  A  and  B  should  be  diagonal  matrices  with  elements  equal  to 
the  inverse  of  the  square  of  the  maximum  allowable  state  excursions  for 
A,  and  the  inverse  of  the  square  of  the  maximum  allowable  control  usage 
for  B.  Although  such  a  procedure  has  been  shown  to  provide  a  reasonable 
first  guess  for  the  matrices,  it  by  no  means  eliminates  the  trial-and- 
error  process.  In  addition,  it  is  really  designed  to  handle  the  steady- 
state  regulator  problem  rather  than  the  meeting  of  transient-response 
sped  f  ications. 

A  different  approach  to  the  problem  is  known  as  destabilization.  In 
this  method,  the  cost  function  is  modified  by  the  multiplication  of  the 


factor  exp(at),  uith  a>0.  Since  the  optimal  control  solution  is  guaran¬ 
teed  to  be  stable,  the  resulting  closed-loop  roots  are  "stabilized"  or 
"shifted  to  the  left"  in  the  s-plane,  by  an  amount  a.  The  problem  uith 
this  approach  is  that  all  of  the  roots  are  shifted,  even  the  ones  uhich 
already  have  an  acceptable  response.  Thus,  control  effort  is  uasted, 
and  some  variables  end  up  being  over-damped  and/or  faster  than  required. 

By  combining  the  techniques  of  model  follouing  and  optimal  control, 
it  is  possible  to  take  advantage  of  the  pouer  and  flexibility  of  the  op¬ 
timal  control  methods  uhile  also  incorporating  time-response  informa¬ 
tion,  resulting  in  a  more  straightforward  design  procedure.  This  proce¬ 
dure,  sometimes  known  as  "Optimal  Model  Follouing",  may  be  accomplished 
using  either  of  two  methods:  state  matching  or  dynamics  matching. 

3.2  STATE  MATCHING 

Let  a  plant  be  described  by 


x  =  Fx  +  Gu 

y  =  Hx, 


uith  the  desired  system  performance  being  given  by  the  model  equations 


X  hi  =  F  nX  m 


y  m  “  Hq|X  m  • 
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The  H  and  Ha  matrices  are  "selector"  matrices  which  identify  which 
states  are  to  be  matched.  Since  we  would  like  to  have  y  be  as  close  as 
possible  to  y m,  we  attempt  to  minimize  the  cost  function 


J  = 


l (ym  ■  y)tA(y.  -  y)  +  u*Bu]dt, 


with  the  system  dynamics  being  given  by 


F  0 
0  F  m 


If  we  define  a  new  state  variable 


then  the  cost  function  may  be  rewritten  as 


J  = 


(z*A*z  +  u*Bu)dt» 
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where 


H  *AH 


-H  *AH* 
Ha^AHn 


The  resulting  control  lau  will  be 


u  =  Cz  =  Cix  +  C2xm, 


and  therefore  requires  both  plant  states  and  model  states  to  implement 
it. 

With  the  state-matching  approach,  it  is  the  "speed"  of  the  plant 
roots  uhich  is  determined  by  the  choice  of  weighting  matrices  in  the 
cost  function.  By  making  these  roots  fast  enough,  the  plant  transfer 
function  approaches  unity.  As  a  result,  the  plant  outputs  will  accu¬ 
rately  "follow"  any  input  signal.  In  this  case,  the  input  signal  is 
simply  the  output  of  our  model  system. 

Unfortunately,  if  the  plant  roots  are  not  made  sufficiently  fast, 
either  because  of  control  limitations  or  for  other  reasons,  they  may  be¬ 
come  the  dominant  factor  in  the  overall  system  response,  overshadowing 
the  mode)  roots.  Thus,  even  if  good  steady-state  performance  is 
achieved,  the  resulting  transient  response  may  be  totally  different  from 
the  desired  one.  In  addition,  since  it  is  not  possible  to  have  the 
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plant  output  exactly  follow  the  input  unless  the  transfer  function  is 

» 

identically  one*  we  would  not  be  able  to  achieve  perfect  model  following 
with  this  type  of  system  unless  the  plant  became  infinitely  fast. 

3.3  DYNAMICS  MATCHING 

Suppose  that  instead  of  trying  to  minimize  the  error  between  the 
plant  and  model  states,  we  attempt  to  match  the  dynamics.  Given  a  plant 

x  =  Fx  +  Gu 

y  =  Hx, 

and  a  model 

we  would  like  the  dynamics  of  y  to 
words,  we  want  y  to  be  as  close  as 
we  can  minimize  the  cost  function 


J  =  jl(y  -  P«y)*A(y  -  F*y)  +  utBuldt, 
0 

which  is  equal  to 


match  the  dynamics  of  xm.  In  other 
possible  to  F„y.  To  accomplish  this 
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(  (hfx  +  HGu  -  FBHx)*A(HFx  +  HGu  -  F,„Hx)  +  u*Buldt. 

0 


Multiplying  out  and  regrouping  terms,  we  get 


J 


{x*[(HF  -  F0H)*A(HF  -  FmH)x  +  2x M CHG ) *A(HF  -  FmH 

. 

0 

+  u*l(HG)tA(HG)  +  Bluldt. 


Defining 


A  =  [(HF  -  F„H)tACHF  -  FmH>) 
S  =  t(HG)*A(HF  -  FmH)l 
B  =  t(HG)*A(HG)  +  Bi, 


ue  can  write 


J 


to 

[(x  lAx  +  2x  *S  *u  +  u*Bu)dt. 


0 
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If  the  cross-product  term  could  be  made  to  disappear,  the  problem  would 
then  be  in  the  standard  optimal  control  form.  Kriendlerf 161  points  out 
that  frequently  the  plant  being  worked  on  will  contain  actuators,  which 
isolate  the  plant  states  of  interest  from  the  control  inputs.  In  such 

A  A 

cases,  HG=0.  Therefore,  S=0  and  B=B,  and  we  are  left  with 

» 

» 

J  =  (x*Ax  +  utBu)dt, 

1 

0 

which  is  easily  solved  in  the  standard  way.  However#  with  a  little  ad¬ 
ditional  manipulation  (as  suggested  by  Anderson  and  Moorelll),  even  this 
restriction  can  be  removed.  Completing  the  square, 


•  A  A  .A  A 

A  =  A  -  S*8*  *S 


and 


A  A 

Ui  -  u  ♦  8* 'Sx. 


It  is  interesting  to  examine  the  case  when  perfect  model  following  is 
possible  (as  discussed  in  Chapter  II),  and  no  weight  is  put  on  the  con¬ 
trol.  For  this  problem  A  becomes  0,  which  in  turn  means  that  no  weight¬ 
ing  will  be  assigned  to  the  states.  As  a  result,  the  output  of  an  opti¬ 
mal  control  program  will  be  2ero  gains  for  Ui  (assuming  the  system  was 
stable  to  start  with).  If  ue  let  A  be  the  identity  matrix  (i.e.,  each 
state  match  is.  equally  important),  then 

A 

B  =  l(HC)*(HG)  1 

A 

S  =  ((HG)*(HF  -  FmH)J, 

and  the  control  is  given  by 

U  =  -l(HG)*(H3)l*,(HG)t(HF  -  F„ll)x 
=  (HG)*(F„H  -  HF)x, 
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the  same  results  which  were  achieved  in  Chapter  II.  although  admittedly, 
this  time  after  a  much  more  complex  procedure.  Thus,  with  the  dynamics¬ 
matching  approach.  ue  can  achieve  perfect  model  following  when  it  is 
possible,  without  the  need  for  an  infinitely  fast  plant. 

3.4  GENERAL  REMARKS 

In  the  preceeding  sections  it  was  shown  how  the  ideas  of  model  fol¬ 
lowing  and  optimal  control  could  be  combined  to  produce  "Optimal  Model 
Following"  design  methods.  With  regular  optimal  control,  selection  of 
the  weighting  matrices  is  a  "black  art"  at  best.  Putting  too  little 
weight  on  the  control  results  in  ridiculously  fast  roots  and  impracti¬ 
cal  ly  high  control  levels,  and  changing  the  weighting  on  one  state  usu¬ 
ally  changes  the  response  of  other  variables  as  well.  With  Optimal 
Model  Following,  the  designer  need  only  formulate  a  model  whose  dynamic 
equations  satisfy  the  time  response  specifications.  Typically,  ue  might 
then  let  A  be  the  identity  matrix  (in  other  words,  all  state  matches  are 
equally  important),  and  let  B  equal  a  very  small  number.  If  perfect 
model  following  is  possible  and  the  dynamics-matching  method  is  used,  B 
can  even  be  set  to  zero.  In  this  case  the  control  gains  will  not  in¬ 
crease  past  the  finite  gains  which  are  required  for  the  match,  no  matter 
what  A  is,  and  the  plant  states  will  exactly  match  the  model  states. 
Kreindlert  16!  noted  that  with  an  actuator  in  the  plant  but  not  in  the 
model,  use  of  a  small  B  led  to  high  feedback  gains  on  the  actuator,  es¬ 
sentially  trying  to  "speed  up"  the  actuator  root.  This  makes  sense  when 
ue  consider  the  perfect-model-fol lowing  ideas  of  Chapter  II.  With  an 
actuator  in  the  system,  some  of  the  plant  states  are  not  directly 
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controllable,  ao  perfect  node!  following  is  not  possible,  in  general. 
Closer  and  closer  matches  can  be  achieved  by  speeding  up  the  response  of 
the  state  through  which  the  plant  states  are  controlled  (i.e.  the  actua¬ 
tor).  However,  because  of  nonlinearities  or  actuator  saturation,  there 
is  usually  a  physical  limit  to  the  actuator  speed  of  response.  In  any 
case,  whether  perfect  model  following  is  possible  or  not,  the  designer's 
task  is  a  very  straightforward  one.  Even  if  some  iteration  is  required, 
it  is  always  obvious  which  direction  to  proceed.  If  the  control  use  is 
too  high,  just  increase  the  B  matrix.  If  the  degree  of  matching  is  too 
low,  just  decrease  B.  In  this  way  the  designer  can  easily  achieve  a 
balance  between  control  use  and  the  "goodness"  of  the  match  of  the  de¬ 
sired  time  response.  The  method  thus  serves  as  a  kind  of  link  between 
the  time-response  ideas  of  classical  control  theory,  and  the  cost  func¬ 
tion  and  "control-oriented"  approach  of  optimal  control  theory. 

What  are  some  of  the  problems  or  shortcomings  of  these  methods?  An 
obvious  disadvantage  of  all  explicit  model-following  systems  is  the  nec¬ 
essity  to  implement  the  model  equations.  This  requirement  means  that 
the  control  computer  will  probably  have  to  be  larger  and  faster  than  it 
otherwise  uould  have  to  be.  On  the  other  hand,  as  we  will  see  in  Chap¬ 
ter  IV,  explicit  model-following  systems  have  advantages  in  disturbance 
rejection  and  control -system  sensitivity.  (State-matching  systems  re¬ 
quire  an  explicit  formulation,  while  dynamics-matching  systems  may  be 
either  explicit  or  implicit.)  If  a  plant  state  is  not  directly  cont¬ 
rollable,  attempts  to  match  the  model  dynamics  closely  will  result  in  a 
speeding  up  of  the  transfer  state  (actuator)  according  to  how  heavy  the 
weighting  on  the  match  was  made.  Finally,  because  the  objective  in 
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optimal  control  problems  is  to  minimize  a  cost  function, 
aluays  be  an  answer,  even  if  it  is  totally  unacceptable.  Thus,  the  use 
of  very  ambitious  models  without  the  willingness  to  expend  large  am¬ 
ounts  of  control,  or  the  attempt  to  match  many  states  with  only  a  feu 
controls,  will  likely  yield  disappointingly  poor  matches.  As  a  result, 
it  is  probably  a  good  idea  to  review  the  perfect-model -fol 1  owing  crite¬ 
ria  to  get  a  feel  for  what  is  possible  before  determining  the  objec¬ 
tives.  And,  once  a  final  design  has  been  arrived  at,  system  testing  or 
simulation  is  a  necessity  for  verifying  the  performance  of  complex  sys¬ 
tems. 

3.5  EXAMPLE :  THO  MASSES  AND  £  SPRING 

As  an  example  of  the  ease  of  the  Optimal  Model  fol lowing  design  pro¬ 
cess,  specifically  the  dynamics-matching  method,  consider  the  system 
i  consisting  of  two  masses  and  a  spring  shown  in  figure  2. 

The  single  control  acts  only  on  the  first  mass,  but  it  is  desired  to 
k  control  the  motion  of  the  second  mass. 

v  The  plant  equations  of  motion  are  given  by 
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Me  would  like  the  second  mass  to  behave  like  a  second-order  system 
with  a  frequency  of  2  rad/sec  and  a  damping  ratio  of  .7.  The  model  is 
therefore  described  by 
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If  ue  were  trying  to  solve  this  problem  with  the  standard  Linear 
Quadratic-  Regulator  (IQR)  approach,  we  might  initially  try  using  an  A 
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Matrix  with  unity  weighting  on  the  position  of  tho  second  mass  and  ze¬ 
roes  elsewhere.  Ue  could  then  try  several  different  values  of  B  to  see 
what  would  happen.  A  plot  of  the  resulting  root  locations  is  given  in 
Figure  3. 

Since  the  cost  function  does  not  contain  any  direct  information  about 
the  desired  transient  response,  it  would  have  been  a  lucky  coincidence 
if  the  plant  roots  had  ended  up  close  to  the  model  roots.  It  is  no 
doubt  possible  to  do  much  better  than  ue  have  done  here,  given  enough 
tries.  Ue  can  increase  the  weighting  on  one  of  the  states  to  speed  up 
that  state,  or  increase  the  weighting  on  the  derivative  of  a  state  to 
increase  its  damping.  However,  this  is  clearly  a  trial-and-error  pro¬ 
cess,  and  it  may  take  a  large  number  of  attempts  before  satisfactory  re¬ 
sults  are  achieved. 

If  ue  use  Optimal  Model  Following,  the  design  process  is  greatly  sim¬ 
plified.  Since  ue  are  trying  to  match  the  third  and  fourth  plant  states 
to  the  model,  the  selector  matrix  is 


H 


C  0  1  0 

0  0  0  1 


The  results  are  shoun  in  Figure  4. 

Lowering  the  weighting  on  the  control  causes  two  of  the  plant  roots 
to  proceed  rapidly  toward  the  model  roots,  while  the  remaining  plant 
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Figure  4:  Locus  of  Roots  v.s.  Control  Weighting  -  OMF  Design 
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roots  become  vory  fast.  Sines  ths  sscond  mass  is  not  directly  controll¬ 
able.  the  first  mass  acts  like  an  actuator  and  becomes  faster  and  faster 
as  the  weighting  on  the  control  is  decreased.  The  designer  needs  only 
to  decide  upon  the  maximum  acceptable  control  use.  which  will  then  det¬ 
ermine  how  closely  the  plant  states  will  come  to  achieving  the  desired 
time  response.  It  should  be  noted  that  Optimal  Model  Following  is  not  a 
pole-matching  technique,  and  therefore  the  plant  roots  may  not  always 
proceed  directly  to  the  model  roots.  However,  the  time  responses  of  the 
plant  states  will  get  closer  and  closer  to  the  time  responses  of  the 
model  as  the  control  weighting  is  decreased.  A  plot  of  the  impulse  res¬ 
ponse  of  the  second  mass  with  an  LQR  design  is  given  in  Figure  5.  while 
the  results  for  two  OMF  designs  are  presented  in  Figure  6. 
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Figure  6«  Impulse  Response  for  OMF  Designs 


Chapter  IV 


tlODEL  FOLLOWING  FOR  INSENSITIVE  CONTROL 

4.1  REVIEW  OF  PAST  EFFORTS 

Unknown  or  varying  parameters  have  always  been  a  concern  to  control- 
system  designers.  Classically,  this  was  handled  by  conservatively  se¬ 
lecting  the  control -system  gains  and  structure  so  that  the  predicted 
gain  and  phase  margins,  damping  ratios,  and  time  constants  indicated 
that  adequate  performance  could  be  achieved  even  for  "worst  case"  condi¬ 
tions.  This  was  then  confirmed  through  simulation  or  system  test,  and 
if  necessary,  changes  were  made  to  the  design. 

On  the  surface,  optimal  control  techniques  appeared  to  make  the  de¬ 
sign  process  much  more  objective.  It  was  only  neoessary  to  input  a  feu 
weighting  parameters  to  the  computer,  and  the  "optimum"  solution  would 
be  calculated.  The  key  question  of  course,  is  "Optimum  with  respect  to 
what  criteria?"  Kalman! 1 4 J  has  shown  that  almost  any  control  law  can  be 
shown  to  be  an  optimum  one,  given  the  proper  choice  of  cost  function. 

The  sensitivity  of  a  control  system  may  be  defined  in  many  different 
ways,  such  as  trajectory  sensitivity,  eigenvalue  sensitivity,  or  cost 
function  sensitivity.  Typically  though,  the  problem  is  that  while  con¬ 
trol  systems  designed  using  standard  linear  quadratic  regulator  tech¬ 
niques  may  perform  very  well  for  the  "nominal"  conditions,  with  only 
small  errors  in  the  parameters,  system  performance  can  be  significantly 


degraded  or  even  unstable.  This  has  no  doubt  led  to  much  confusion  and 
embarrassment  among  newly  graduated  engineers.  trying  to  apply  the 
theory  they  learned  in  school  to  real-world  problems.  It  has  probably 
also  helped  to  prevent  more  widespread  use  of  modern  control  techniques 
in  the  aerospace  industry.  Actually,  sensitivity  is  not  an  inherent 
characteristic  of  the  optimal  control  approach,  but  rather  a  result  of 
the  designer  failing  to  incorporate  all  of  his  design  requirements  into 
the  criteria.  As  we  shall  see  in  a  later  example,  the  standard  LQR 
theory  can  usually  be  used  very  successfully  to  give  designs  with  as 
much  stability  as  desired,  even  for  large  parameter  uncertainties,  as 
long  as  the  designer  takes  these  uncertainties  into  account  during  the 
design  process.  Although  the  resulting  system  no  longer  behaves  "opti¬ 
mally"  at  the  nominal  parameter  values,  at  least  the  instability  problem 
is  avoided. 

Over  the  past  few  years,  there  have  been  a  number  of  attempts  to 
somehow  automatically  incorporate  sensitivity  characteristics  into  the 
design  process.  Palsson  and  Uhitaker[23]  and  Peled[24]  applied  random- 
vector  approaches  with  some  success.  Harvey  and  PopeI13J  and 
Vinkler[30]  produced  excellent  comparisons  of  a  number  of  the  most  com¬ 
mon  methods  for  sensitivity  reduction.  Although  some  of  the  methods  ap¬ 
pear  to  be  better  than  others,  no  single  technique  has  been  found  to  be 
clearly  superior  to  the  others  in  all  cases.  Almost  all  of  the  methods 
shared  two  common  characteristics.  First,  computational  requirements 
were  very  great.  Most  of  the  algorithms  required  the  use  of  gradient- 
search  routines  to  find  the  minimum  of  a  cost  function,  and  many  itera¬ 
tions  were  sometimes  required.  Second,  and  most  importantly,  system 
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performance  and/or  control  ute  at  the  nominal  parameter  values  was 
degraded,  with  the  amount  depending  on  hou  insensitive  the  system  uas 
asked  to  be.  In  addition,  although  stability  uas  guaranteed  in  some  of 
the  methods,  system  eigenvalues  could  vary  dramatically  with  changes  in 
the  parameters.  The  model-fol lowing  method  discussed  later  in  this 
chapter  is  shown  to  avoid  these  undesirable  characteristics. 

4.2  FEEQFQRUAR D/FEEDBACK  EQUIVALENCE 

Suppose  we  have  a  plant 


*  =  Fx  +  Gu, 


and  a  model,  of  the  same  dimension. 


F»xm  ♦  Gnum. 


Let  us  implement  an  explicit  model-following  control  system,  with  con 
trol  lau  of  the  form 

u  =  C,x  +  Cjx*  +  C3um. 


The  system  structure  is  then  as  shown  in  Figure  7. 

By  adding  and  subtracting  Cjx,  the  control  law  may  be  rewritten  as 
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Figure  7:  Explicit  Model  Following  System 


u  =  CiX  +  C2xn  +  (C2x  -  Cjx)  +  C3uB. 


Thus, 


u  =  (Ci  +  C2)x  +  C2(xm  -  x)  +  C3u„. 


Defining  »=xm-x,  we  obtain 


u  =  (Ci  +  C2)x  +  C2e  +  C3um. 
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The  block  diagram  in  Figure  8  can  thus  be  seen  to  be  equivalent  to  the 
one  in  Figure  7. 


Suppose  the  model  and  the  closed-loop  plant  are  exactly  the  same,  and 
the  initial  conditions  are  also  identical.  Then  x  and  xm  will  aluays  be 
equal,  and  the  error,  e,  will  aluays  be  zero.  In  this  case,  as  Figure  8 
makes  clear,  it  is  (Ci+Cz)  which  has  determined  the  closed-loop  plant 
dynamics.  The  match  of  plant  and  model  states  could  therefore  have  been 
achieved  uith  only  feedback  from  the  plant  states  (Cj=0),  with  only 
feedforward  from  the  model  states  t C i =  0 ) ,  or  uith  any  combination  of  the 
two.  as  long  as  the  sum  of  Ci  and  Cz  remained  the  sane. 
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As  the  above  discussion  indicates*  Implicit  Model  Following  is  just 
as  effective  as  Explicit  Model  Following  in  terms  of  matching  dynamics 
(assuming  plant  and  model  are  of  the  same  dimension  and  all  of  the 
states  are  matched).  Thus,  in  cases  where  the  plant  parameters  are  well 
known  and  there  are  no  unknown  disturbances,  or  where  a  conservative  de¬ 
sign  is  acceptable,  there  is  no  need  to  implement  the  model.  Where  this 
is  not  the  case,  the  explicit  formulation  should  be  used,  as  discussed 
in  the  following  sections. 

4.3  DESCRIPTION  £F  IH£  DESIGN  METHOD 

Many  authors  have  stated  their  "intuitive"  feelings  that  explicit 
model-following  systems  should  have  some  advantages  in  terms  of  system 
sensitivity,  at  least  as  compared  to  Implicit  Model  Following  systems. 
It  certainly  seems  reasonable  that  the  ability  to  generate  an  error 
signal  between  the  plant  and  a  model  system  and  using  it  as  part  of  the 
control  input  would  help  to  keep  plant  states  "close"  to  the  model 
states.  However,  unless  the  system  roots  are  selected  carefully,  the 
transient  response  of  the  plant  may  bear  no  resemblence  to  the  model 
response . 

Tylerl29]  and  others  formulated  the  explicit  model-following  problem 
as  finding  the  control  which  minimized  the  cost  function 


J  = 


(e*Ae  +  u*Bu)dt. 
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Thus,  they  reasoned  that  an  optimal  tradeoff  could  be  accomplished  bet¬ 
ween  the  error  and  the  control  effort  used. 

Actually,  the  above  cost  function  uas  usually  computed  using 


J 


x ) *A ( x „  -  x)  +  u*Buldt. 


As  a  result,  optimal  gains  were  computed  for  both  x  and  xm.  even  though, 
as  ue  sau  in  the  proceeding  section,  the  model-fol lowing  ability  of  a 
system  is  not  dependent  on  whether  it  is  implemented  with  feedforward, 
feedback,  or  a  combination  of  the  two.  In  other  words,  the  gains  which 
minimize  the  above  cost  function  do  not  uniquely  determine  the  plant 
response,  and  different  combinations  of  feedforward  and  feedback  gains 
will  have  substantially  different  effects  for  off-nominal  conditions. 

• 

Kriendlerl 161  noted  the  poor  dynamics-matching  abilities  of  this  t/pe 
of  system,  uhich  caused  him  to  prefer  Implicit  Model  Following.  The 
poor  matches  really  had  nothing  to  do  with  whether  the  system  was  expli¬ 
cit  or  implicit,  but  rather  was  a  function  of  using  state  matching  in¬ 
stead  of  dynamics  matching.  Kriendler  also  acknowledged  that  Explicit 
Model  Following  might  have  better  sensitivity  properties,  although, 
since  he  used  Tyler's  method  (state  matching)  to  calculate  the  gains,  he 
found  no  significant  difference  in  sensitivity  between  the  two  methods. 
In  ar.y  event,  to  incorporate  the  advantages  of  both  explicit  and  impli¬ 
cit  model-following  systems,  he  proposed  (although  he  did  not  inveri’- 
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gate)  the  cost  function 


J 


I(y  -  F»y)tA1(y  -  F,y)  ♦  (y„  -  y^A^y,  -  Vi  *  u*Bu]dt 

0 

This  does  appear  to  represent  a  great  improvement  over  the  standard 
state-matching  formulation.  The  error-rate  term  tends  to  match  the  dy¬ 
namics,  while  the  error  term  tends  to  speed  up  the  plant  states  by  in¬ 
creasing  the  feedback  gains.  This  increased  feedback  tends  to  desensi¬ 
tize  the  system  to  parameter  errors.  The  disadvantage  is  that  changes 
in  A2  to  adjust  system  sensitivity  will  also  result  in  changes  to  the 
nominal  plant  dynamics/transient  response.  Instead,  it  is  claimed  that 
the  design  process  should  be  accomplished  twice:  once  to  establish  the 
plant  dynamics  at  the  nominal  parameter  values,  and  once  to  establish 
the  error  dynamics  (and  therefore  the  system  sensitivity).  By  doing  the 
design  in  this  manner,  system  performance  at  the  nominal  parameter  val¬ 
ues  becomes  independent  of  the  sensitivity. 

The  proposed  design  method,  which  we  call  Model  Following  for  Insen¬ 
sitive  Control,  is  then  as  follows: 

1.  Perform  a  design,  using  any  method  (root  locus,  model  follow¬ 
ing,  optimal  control,  etc.),  to  select  the  nominal  plant  res¬ 
ponse.  These  gains  become  Ct=Ci+Ci  and  C3. 

2.  Implement  a  model  of  this  closed-loop  system.- 

3.  Oo  another  design  to  select  the  error  roots,  which  should  be 


chosen  faster  than  the  plant  roots.  These  gains  become  Ci. 


4.  Calculate  CZ=C t~Ci . 

5.  Implement  the  control  law  as  u=Cix+C2Xb+C3U«. 

4.4  THEORETICAL  DEVELOPMENT 

Let  the  behavior  of  a  plant  with  nominal  parameter  values  be  given  by 

x  =  Fx  +  Gu. 

However,  due  to  parameter  uncertainty  and/or  slowly  varying  parameters, 
suppose  the  plant  actually  responds  according  to 

x  =  CF  +  AF)x  +  Gu. 

Through  the  use  of  model -fol lowing  techniques,  optimal  control,  root  lo¬ 
cus,  or  any  other  method,  assume  that  a  satisfactory  closed-loop  control 
lau  can  be  obtained  for  nominal  values  of  the  parameters.  Let  this  con¬ 
trol  Ir.u  be  given  by 

u  =  C  tx  +  Cju*. 

Then,  for  nominal  values  of  the  parameters,  we  have 

x  =  Fx  +  G(C*x  +  C3Uffl) 
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=  (F  +  ec*)x  ♦  GC3u«. 


Build  a  math  model  of  this  system,  i.e., 

x,,  =  (F  +  GCt)xm  +  GC3u« 

=  FpXji  ♦  GnUin. 

For  the  second  design,  ue  again  may  use  whatever  method  we  desire,  but 
the  objective  this  time  is  to  select  the  error  roots.  For  an  insensi¬ 
tive  design,  these  roots  should  be  chosen  faster  than  the  nominal-plant 
roots.  The  resulting  gains  are  designated  Ci . 

I4e  then  use  a  control  law 

u  =  C,x  +  C2xm  +  C3um 

=  C*x  +  C2e  +  C5Um, 

where  e=xm-x  and  Ct=Ci+C2. 

This  results  in 

x  =  (F  +  AHx  +  G(Ctx  +  C2e  +  C3u  m) 

=  (F  +  AF  +  GCt)x  +  CC2e  +  GC3um. 
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Subtracting,  ue  gat 


e  =  CF  +  GCt)e  -  AFx  -  GC2e 
=  (F  +  GC \  -  GCZ )e  -  AFx . 


The  system  is  thus  described  by 


X 

F+AF+GC  t  GC2 

X 

gc3 

r 

+ 

e 

-AF  F+GC  t“GC2 

e 

0 

uith  a  characteristic  equation  given  by 


sI-F-AF-GC*  -GC2 
AF  si  -F-GC  v*-GC2 


0. 


Since  the  determinant  of  a  matrix  is  not  changed  by  elementary  rou  and 
column  operations,  ue  can  urite 


Sl-F-GCt  sI-F-GCt 
AT  sI-F-GC  t+GC* 


0 


and 
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Sl-F-GCt 


0 


=  0. 


jar  si-r-AF-BCt*GC2l 


This  can  be  written  as 


(Si  -  F  -  GCtXsI  -  F  -  AF  -  GC*  +  CC2 >  =  0, 


or  equivalently. 


(si  -  Fm)(sl  -  F  -  AF  -  GC,)  =  0. 


Thus,  given  an  appropriate  choice  for  <o,  the  roots 
can  usually  be  made  fast  enough  to  keep  the  roots  of  (si 
inant  system  roots,  even  for  arbitrarily  large  AF.  If 
are  made  sufficiently  fast,  they  will  not  significantly 
tern  response,  regardless  of  the  value  of  AF.  Also, 
model  start  with  identical  initial  conditions  and  AF=0, 
will  never  even  be  excited,  since  the  system  equations  w 


X 

F+GC  t 

G  C  2 

X 

* 

GC  j 

z 

+ 

• 

e 

0 

F+GC Z 

e 

0 

f  (sI-F-AF-GC,) 
Fm)  as  the  dum- 
the  error  roots 
affect  the  sys- 
f  the  plant  and 
the  error  roots 
11  be 


4.5  EXAMPLE;  FIRST-ORDER  SYSTEM 


In  order  to  get  a  better  understanding  of  how  Model  Following  for  in¬ 
sensitive  Control  works,  let  us  consider  the  simple,  first-order  system 
used  by  Shenkar(26I.  The  plant  is  described  by 


x  =  - ( 1  +  AF)x  +  u. 


where  AF  represents  the  uncertainty  in  a  plant  parameter.  The  value  of 
AF  may  be  anything  from  -1  to  1,  with  0  being  the  nominal  value.  Three 
different  design  techniques  were  evaluated:  the  standard  Linear  Quad¬ 
ratic  Regulator  (LQR)  technique,  the  Average  Cost  Function  ( ACF )  ap¬ 
proach  (representative  of  the  methods  of  Vinkler,  Ashkenazi,  and 
Shenkar),  and  Model  Following  for  Insensitive  Control  (MriC). 

For  the  LQR  design,  we  assume  that  our  design  requirements  are  met  by 
letting  A=B=1  for  the  nominal  plant.  However,  as  might  be  expected,  the 
system  dynamics  are  very  dependent  on  the  value  of  the  uncertain  parame¬ 
ter.  The  impulse  responses  for  the  nominal  and  extreme  values  of  AF  are 
shown  in  Figure  9. 

With  the  Average  Cost  Function  design  there  is  definitely  some  im¬ 
provement,  but  we  still  have  a  wide  range  of  responses  for  the  different 
values  of  AT.  Also,  even  if  AF=0,  the  response  is  now  different  from 
the  nominal  LQR  design.  Results  for  this  design  are  shoun  in  Figure  10. 
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-Order  System  -  IQR  Design 


-Order  System  —  ACF  Dtsifln 


For  the  Model  Following  for  Insensitive  Control  design,  the  closed- 
loop  root  for  the  nominal  plant  was  chosen  to  be  the  same  as  the  nominal 
LQR  system,  while  the  error  root  was  arbitrarily  chosen  to  be  at  s=-10. 
With  this  system,  the  impulse  responses  are  almost  identical,  regardless 
of  the  value  of  AF .  The  responses  are  shown  in  Figure  11. 

If  a  more  insensitive  design  is  required,  one  can  simply  select  a 
faster  error  root.  Without  disturbances,  the  control  levels  required 
are  directly  related  to  the  transient  response.  Thus,  the  control  mag¬ 
nitudes  for  the  MF I C  design  are  essentially  the  same  as  the  control  mag¬ 
nitudes  which  would  have  been  required  for  a  design  in  which  AF  was 
known  in  advance.  Control  bandwidth  however,  obviously  must  be  in¬ 
creased  as  the  error  roots  are  speeded  up  (which  may  be  costly). 

To  understand  what  is  happening  physically,  it  is  helpful  to  think  in 
terms  of  the  system  transfer  functions.  The  open-loop  plant  is  des¬ 
cribed  by  1/Cs+l).  The  desired  closed-loop  system,  which  is  attained 
with  the  nominal  LQR  design,  has  a  transfer  function  of  1/(s+ 1 . 41 4 ) . 
The  ACF  design  has  a  slightly  higher  feedback  gain,  resulting  in  a 
transfer  function  of  1/(s+l.G89).  With  the  MFIC  design,  the  plant  is 
greatly  speeded  up,  in  this  case  to  l/Cs+lO).  However,  the  feedforward 
from  the  model  and  the  model  input  form  a  prefilter  or  "command  condi¬ 
tioner"  with  a  transfer  function  of  Cs+ 1 0 )/(s+ 1 . 4 1 4) .  Thus,  with  nomi¬ 
nal  parameters,  the  speeded-up  plant  root  is  exactly  cancelled,  result¬ 
ing  in  the  desired  system  transfer  function.  With  a  non-zero  AT, 
perfect  cancellation  is  not  attained,  but  because  the  pole  and  zero  are 
close  to  each  other  and  both  are  much  foster  than  the  rest  of  the  sys¬ 
tem,  they  do  not  significantly  affect  the  response. 
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The  beneficial  effects  of  this  design  method  stem  from  tuo  factors: 
first*  the  increased  feedback  gains*  uhich  have  a  desensitizing  effect 
uell-knoun  to  classical  designers*  and  second*  the  cancellation  of  the 
speeded-up  plant  roots.  Note  that  either  of  these  items  by  itself  uould 
not  yield  the  desired  results.  Just  increasing  the  feedback  gains  uould 
louer  the  system  sensitivity*  but  the  system  response  uould  be  over¬ 
damped  or  faster  than  desired.  On  the  other  hand,  just  cancelling  the 
plant  roots  uould  allou  us  to  achieve  the  desired  closed-loop  roots,  but 
it  does  nothing  for  the  sensitivity  problem. 

The  fine  design  procedure  is  primarily  oriented  to  obtaining  a  desir¬ 
able  transient  response,  rather  than  to  disturbance  rejection.  On  the 
other  hand*  disturbance  rejection  is  a  very  important  part  of  the  Aver¬ 
age  Cost  Function  approaches.  To  try  to  compare  the  tuo  methods  from  a 
different  perspective,  the  same  first-order  example  uas  used,  but  uith 
the  addition  of  a  process-noise  source,  uith  a  pouer-spectral-densi ty 
and  a  process-noise  distribution  element  of  1.  The  tuo  components  of 
the  cost  function  -  state  excursions  and  control  use  -  uere  then  plotted 
against  feedback  gain  in  Figures  12  and  13. 

As  might  be  expected,  the  greater  the  negative  feedback  gain,  the 
tighter  the  control  system  is  able  to  regulate  the  state,  and  hence  the 
smaller  the  state  excursions.  Likewise,  the  greater  the  feedback  gains 
are.  the  higher  the  control  level  that  is  used  to  handle  the  distur¬ 
bance.  By  adding  these  tuo  terms  together,  the  total  value  of  the  cost 
function  is  obtained.  Plots  of  this  total  cost  function  versus  the  un¬ 
certain  parameter  are  presented  for  the  LQR  design,  the  ACF  design,  and 
tuo  HFIC  designs  in  Figure  14. 
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The  advantage  of  the  ACF  design  is  supposed  to  be  that  by  accepting  a 
slightly  higher  value  of  the  cost  function  for  positive  and  zero  values 
of  the  uncertain  parameter,  a  substantial  improvement  can  be  achieved 
for  negative  values  of  the  parameter.  However,  the  MFIC  design,  which 
had  much  better  transient  responses,  shows  significantly  higher  values 
of  the  cost  function  (which  contains  control  effort)  for  any  &F.  in 
fact,  the  faster  the  error  root  is  made,  (and  hence  the  more  insensitive 
the  resulting  system),  the  higher  the  cost  function  becomes. 

The  message  from  this  example  is  that  it  is  important  for  the  de¬ 
signer  to  be  sure  that  he  has  chosen  the  cost  function  which  properly 
represents  the  problem  he  is  trying  to  solve.  In  this  case,  by  using  a 
model-following  structure  and  increased  feedback  gain,  it  was  possible 
to  achieve  a  significant  improvement  in  transient  response  compared  to 
the  so-called  "optimal"  design.  The  cost  of  this  improvement  was  the 
increased  control  usage  due  to  the  process  noise,  which  may  or  may  not 
be  important  to  the  designer, 

4-6  EXAMPLE:  SECOND-ORDER  SYSTEM 


For  our  next  example,  consider  the  second-order  system  studied  by 


F 


0  1 
-2+u^  I+02 


6  = 


0 

1 


and 


-1  i  w,  i  1 
-1  i  ui  i  1.5. 

With  weighting  matrices  of 


A 


1  0 
0  1 


and  B  =  10, 


a  standard  LQR  design  yields  an  unstable  system  for  some  values  of  w. 

Vinkler  evaluated  several  design  techniques  which  have  been  suggested 
for  handling  uncertain  parameters.  He  then  described  tuo  methods  for 
obtaining  better  results:  the  Multistep  Guaranteed  Cost  Control  method 
(MGCC )  and  the  Modified  Discrete  Expected  Cost  method  (MDEC).  Both 
methods  require  substantial  computational  effort  and  special  computer 
programs  during  the  design  stage.  Actually,  in  many  problems  the  de¬ 
signer  can  obtain  satisfactory  results  with  standard  LQR  techniques,  ns 
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long  as  they  art  used  property.  If  the  designer  is  aware  that  the 
system  has  uncertain  parameters,  he  should  make  his  nominal  design  more 
stable  than  he  ordinarily  uould.  Then,  even  in  the  presence  of  destabi¬ 
lizing  values  of  the  parameters,  an  adequate  degree  of  damping  is  still 
obtained.  For  example,  in  the  present  problem,  if  we  choose 


A 


1  0 
0  50 


and  B  =  10. 


the  closed-loop  response  is  stable  for  all  allowable  values  of  u.  The 
disadvantage  with  this  approach  (as  well  as  with  Vinkler's)  is  that  in 
order  to  get  a  satisfactory  stability  level  with  "unfavorable"  values  of 
the  uncertain  parameter,  we  have  to  accept  overdamped  or  faster  than  de¬ 
sired  responses  for  nominal  or  "favorable"  values  of  the  uncertain  par¬ 
ameter.  By  using  Model  Following  for  Insensitive  Control,  we  can  avoid 
this  result.  For  this  example,  the  nominal-plant  roots  were  chosen  to 
be  the  same  as  for  the  original  IQR  design.  The  error  roots  were  arbi¬ 
trarily  chosen  to  both  be  at  s=-5.  The  resulting  eigenvalues  for  each 
method  are  shown  in  Table  1. 

Notice  that  with  the  MFIC  design  the  dominant  eigenvalues  remain  con¬ 
stant,  and  in  fact,  are  identical  to  the  desired  nominal  IQR  roots.  If 
a  more  insensitive  design  is  desired,  faster  error  roots  can  be  selected 
(with  correspondingly  increased  control  effort),  so  that  their  effect  on 


the  system  response  becomes  negligibly  small. 


4.7  EXAMPLE  =  SIMULATION  OF  A  SUPERSONIC  TAMSEQM 


As  s  final  example,  we  apply  the  fIFIC  design  method  to  the  problem  of 
simulating  the  longitudinal  motions  of  a  supersonic  transport.  This 
problem  has  appeared  at  least  twice  in  the  technical  literature:  Uinsor 
and  Roy l 31 ]  used  it  to  demonstrate  .a  trajectory  sensitivity  method,  and 
Landau  and  CourtiolUS!  gave  a  solution  for  an  adaptive-control  ap¬ 
proach.  The  plant  aircraft  is  the  Total  In-Flight  Simulator  (TIFS).  a 
modified  Convair  C-131B  which  has  side-force  generators  and  direct-lift 
flaps  in  addition  to  the  normal  control  surfaces. 

The  plant  equations  of  motion  are  given  by 


x  =  Fx  +  Gu, 

with 

x*  =  19  q  «  v] 


0. 

1.0 

0. 

0. 

1.401x10-" 

f1q 

-1.9513 

.0133 

-2.505x10’" 

1.0 

-1.3239 

-.0235 

-.5610 

0. 

.  35S0 

-.0279 

0. 

0. 

0. 

-5.3307 

6. 447x10° 

-.2669 

-.1600 

-1.155x10** 

-.2511 

0  . 

.1060 

.0862 
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and 


u*  =  tSe  St  Sf ] . 

The  state  vector  consists  of  pitch  angle,  pitch  rate,  angle  of  attack, 
and  velocity.  The  control  vector  is  made  up  of  elevator  deflection, 
throttle  position,  and  flap  deflection. 

The  model  equations  are 

X„  =  F  niX  n  +  6  mU  ib  , 

ui  th 


X n,  =  1 0 in  dm 


0. 

1.0 

0. 

0. 

5. 31 8x  10' 7 

-.4179 

-.1202 

2. 319x10' 3 

-4.619x10-7 

1.0 

-.7523 

-2.387x10'* 

-.5614 

0. 

.  3002 

-1.743x10'* 

0. 

0. 

-. 1717 

7.451x10'* 

-.0238 

-7.783x10'* 

0. 

3. 655x 10' 3 

and 
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U«*  =  is««  fit.]. 


The  parameter  Mq  has  a  nominal  value  of  -2.038,  but  it  is  assumed  to 
vary  from  -.558  to  -3.558,  which  represents  a  variation  of  approximate! y 
75%  around  the  nominal  value.  Since  perfect  model  following  is  possible 
for  this  system,  the  plant  will  exactly  match  the  model  responses  for 
nominal  Mq.  However,  with  Mq=-3.558,  the  responses  change  signifi¬ 
cantly.  By  using  the  MFIC  design  method,  we  can  obtain  much  better  re¬ 
sults.  Time  histories  for  a  step  elevator  input  for  the  model,  the 
off-nominal  plant,  and  the  off-nominal  plant  using  a  MFIC  design  are 
given  in  Figures  15-17.  The  error  roots  were  arbitrarily  chosen  to  bo 
at  s=-10,  -10,  -2,  and  -2,  with  the  faster  roots  corresponding  to  the 
pitch-angle  and  pitch-rate  errors.  As  shown  in  Figure  18,  the  maximum 
pitch-angle  error  was  less  than  .005  degrees.  The  angl e-of-attack  and 
velocity  responses  had  zero  error.  This  represents  a  considerable  im¬ 
provement  over  the  tra jectory-sensi ti vi ty  method  tried  by  Winsor  and  Roy 
and  the  adaptive  design  given  by  Landau  and  Courtiol .  Control  uses  are 
reasonable  and  are  presented  in  Figures  19-21. 

The  system  response  to  an  initial  error  between  the  plant  and  the 
model  pitch  angle  was  also  investigated,  with  nominal  values  for  the 
parameters.  The  results  are  shown  in  Figure  22.  The  error  is  driven 
rapidly  to  zero,  and  unlike  the  previously  published  efforts,  errors  are 
not  excited  in  the  angle  of  attack  and  velocity  states. 
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figure  17:  C-131/SST  Velocity  Response 
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Figure  19:  C-131/SST  Elevator  Deflection  with  MFIC 
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figure  20t  C-131/SST  Map  Deflection  vith  MFIC 


Figure  21 i  C-131/SST  Throttle  Position  with  nric 


«•*  LIMITATIONS 


The  Model  Following  for  Insensitive  Control  design  method  is  based  on 
the  beneficial  effects  of  increased  feedback  gains  on  system  sensitiv¬ 
ity.  As  a  result/  it  is  only  practical  for  those  systems  which  have 
sufficient  control  authority  and  control  bandwidth  to  allow  the  plant  to 
be  made  "well-damped"  and  "substantially  faster"  than  the  desired  final 
system  characteristics. 

The  problems  studied  here  have  all  assumed  full-state  feedback.  As 
has  been  shown,  many  such  systems  have  quite  different  responses  and 
even  instabilities  in  the  presence  of  parameter  uncertainties.  Another 
large  class  of  control-system  sensitivity  problems  has  to  do  with  the 
estimator  design.  Some  work  has  been  done  in  this  area  (see  for  example 
(91),  with  the  basic  conclusion  being  that  greater  insensitivity  can  be 
achieved  by  speeding  up  some  of  the  estimator  roots.  The  penalty  for 
doing  this  is  a  decrease  in  the  filtering  action  of  the  estimator,  thus 
resulting  in  a  noisier  output. 

Only  uncertainties  in  the  dynamics  matrix  were  addressed  in  this  re¬ 
search.  No  attempt  has  been  made  to  desensitize  systems  to  errors  in 
the  control  matrix,  and  additional  work  is  needed  in  this  area. 


Chapter  V 


CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 

5.1  CONCLUSIONS 

Model  following  can  be  an  extremely  useful  tool  for  the  control -sys¬ 
tem  designer.  The  ideas  of  perfect  model  following  can  provide  simple, 
algebraically  determined  control  laws  for  some  problems  and  a  useful  un¬ 
derstanding  of  system  structure  for  almost  any  application.  By  combin¬ 
ing  model  following  with  optimal  control  theory,  the  process  of  select¬ 
ing  weighting  matrices  is  greatly  simplified  and  a  very  straightforward 
design  procedure  results.  Model  following  can  also  be  used  to  design 
control  systems  which  provide  protection  against  the  adverse  affects  of 
parameter  uncertainties.  The  Model  Following  for  Insensitive  Control 
design  method  allows  the  designer  to  select  the  insensitivity  and  dis¬ 
turbance-rejection  characteristics  of  the  design  independently  of  the 
no-disturbance,  nominal -parameter  performance.  It  takes  advantage  of 
the  favorable  effects  of  high  feedback  gains,  and  then  cancels  the 
speeded-up  plant  roots  by  using  a  model  system  as  a  prefilter.  The 
method  is  most  effective  for  systems  in  which  the  plant  can  be  made  to 
be  "well-damped"  and  "substantially  faster"  than  the  responses  required 
by  the  model  to  be  followed. 
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5.2 


IflE  FUTURE  research 


1.  The  parameter  uncertainty  problems  inveatigated  here  have  as¬ 
sumed  full-state  feedback;  however*  many  controt  system  sensi¬ 
tivity  problems  are  related  to  the  design  of  the  estimator, 
nodel  following,  and  in  particular  the  HFIC  design  method, 
should  be  studied  for  possible  application  to  estimator,  and 
control ler-plus-estimator  problems. 

2.  Additional  research  is  needed  on  the  best  way  to  handle  uncer¬ 
tainties  in  the  control  matrix,  which  is  a  fundamentally  dif¬ 
ferent  problem  than  the  one  addressed  here. 

3.  The  problems  examined  in  this  research  involved  errors  in  the 
dynamics  matrix  of  a  system  with  a  known  structure.  node! 
following  may  also  provide  benefits  in  problems  where  higher- 
order  structural  modes  have  been  ignored,  or  are  not  well 
known.  Further  research  should  be  done  to  investigate  this 
possibility. 


87 


REFERENCES 


1.  Anderson,  B.O.O.,  and  Moore,  J.B.,  Linear  Optimal  Control ,  Prentice 

Halt,  Engleuood  Cliffs,  New  Jersey,  1971. 

2.  Ashkenazi,  A.,  and  Bryson,  A. E.  Jr.,  "The  Synthesis  of  Control 

Logic  for  Parameter  Insensitivity  and  Disturbance  Attenuation," 
paper  S0-1708,  AIAA  Guidance  and  Control  Conference,  Danvers, 
Massachussetts,  August,  1980. 

3.  Asseo,  S.J.,  "Application  of  Optimal  Control  to  Perfect  Model 

following."  Journal  of  Aircraft.  Vol .  7,  No.  4,  July  1970, 
p. 308-313. 

4.  Bryson,  A.E.  Jr.,  and  Ho,  Y.C.,  Applied  Optimal  Control .  Ginn  and 

Company,  Ualtham,  Massachusetts,  1969. 

5.  Chan,  Y.T.,  "Perfect  Model  Following  with  a  Real  Model,"  Preprints 

of  the  Joint  Automatic  Control  Conference,  Columbus,  Ohio,  June 
1973.  p.287-293. 

6. *  Cliff,  E.M.  and  Lutze,  F.H.,  "Application  of  Geometric  Decoupling 

Theory  to  Synthesis  of  Aircraft  Lateral  Control  Systems," 

Journal  of  Aircraft.  Vol.  9,  No.  11,  November  1972,  p.770-776. 

7.  Cliff.  E.M.,  and  Lutze,  F.H.,  "Reply  to  Comment  by 'Paul  R.  Motyka 

and  G.  Warren  Hall,  Journal  of  Aircraft.  Vol.  10,  No.  11, 
November  1973,  p.701. 

8.  Curran,  R.T.,  "Equicontrol 1 abi 1 i ty  and  the  Model -Fol 1  owing 

Problem,"  Ph.O.  Dissertation,  Stanford  University,  July  1971. 

9.  Doyle,  J.C.,  and  Stein,  G.,  "Robustness  with  Observers,"  IEEE 

Transactions  on  Automatic  Contro 1  .  Vol.  AC-24,  No.  4,  August 
1979,  p.607-61  1  . 

10.  Erzberger,  H.»  "Analysis  and  Design  of  Model  Following  Control 

Systems  by  State  Space  Techniques,"  Preprints  of  the  Joint 
Automatic  Control  Conference.  Ann  Arbor,  Michigan,  June  196C, 
p . 572-531 . 

11.  Gran,  R.,  Berman,  H.,  and  Rossi,  M.,  "Optimal  Digital  Flight 

Control  for  Advanced  Fighter  Aircraft,"  Journal  of  Aircraft . 

Vol.  14,  No.  1,  January  1977,  p.32-37. 


88 


12.  Hadass.  Z.«  and  Powell,  J.O..  "Design  Method  for  Minimizing 

Sensitivity  to  Plant  Parameter  Variations. "  AIAA  Journal.  Vot. 

13.  No.  10.  October  1975.  p.1295-1303. 

13.  Harvey.  C.A..  and  Pope.  R.E.,  "Study  of  Synthesis  Techniques  for 

Insensitive  Aircraft  Control  Systems."  NASA  CR-2803.  April  1977. 

14.  Kalman.  R.E.,  "When  Is  a  Linear  Control  System  Optimal?".  Journal 

of  Baa  1c  Engineering  (Transactions  of  the  ASHE).  March  1964. 
p. 51-60. 

15.  Kamlya.  Y.,  "Construction  of  a  Lou  Parameter  and  Disturbance 

Sensitivity  System  by  s  Model  Following  Method."  International 
Journal  of  Control .  Vol .  23.  No.  4,  April  1976,  p.515-524. 

16.  Kreindler.  E.,  and  Rothschild.  D.»  "Model-Following  in  Linear- 

Quadratic  Optimization,"  AIAA  Journal.  Vol.  14,  No.  7,  July 
1976,  p.835-842. 

17.  Kriechbaum,  G.K.L.,  and  Stineman.  R.U.,  "Design  of  Oesirable 

Airplane  Handling  Qualities  via  Optimal  Control,"  Journal  of 
Aircraft.  Vol.  9.  No.  5,  May  1972,  p.365-369. 

18.  Landau.  I.D.,  and  Courtiol.  B.,  "Adaptive  Model  Following  Systems 

for  Flight  Control  and  Simulation,"  Jounal  of  Aircraft  Vol.  9, 
No.  9,  September  1972,  p.668-674. 

19.  Markland,  C.A.,  "Design  of  Optimal  and  Suboptimal  Stability 

Augmentation  Systems."  AIAA  Journal.  Vol.  8,  No.  4,  April  1970, 

P . 673-679. 

20.  Morrow,  L.D.,  and  Bal asubramanian,  R.,  "Real  Model  Following 

Control."  Journal  of  Aircraft.  Vol.  12,  No.  12,  December  1975, 
p . 996-998. 

21.  Motyka,  P.,  and  Hall,  G.,  "Comparison  of  Geometric  and  Response 

Feedback  Approaches  to  Aircraft  Lateral -Di rectional  Decoupling," 
Journal  of  Ai rcraf t .  Vol.  10,  No.  11,  November  1973,  p.69C-7C0. 

22.  Motyka,  P.R.,  Rynaski,  E.G.»  and  Reynotds,  P.A.,  "Theory  and  Flight 

Verification  of  the  TIFS  Model-Following  System,"  Journal  ai 
Aircraft.  Vol.  9,  No.  5,  May  1972,  p.347-353. 

23.  Palsson,  T.,  and  Whitaker,  H.P.,  "Parameter  Uncertainties  in 

Control  System  Design,"  Preprints  ai  lll£  >loint  Automatic  Control 
Conference .  Stanford,  California.  June  1972.  p. 2 -4 f. — 2 5 *4 . 

24.  Peled,  U.,  and  Powell,  J.D.,  "The  Effect  of  Prefilter  Design  on 

Sample  Rate  Selection  in  Digital  Flight  Control  Systems,"  paper 
78-1308,  AIAA  Guidance  and  Control  Conference,  Palo  Alto, 
California,  August  1978. 


89  - 


f 


25.  Rodlsss,  H. A. •  and  Whitaker,  H.P.,  "A  Nau  nodal  Performance  Index 

for  Engineering  Design  of  Flight  Control  Systems, "  Journal  gf 
Aircraft.  Vol .  7,  Ho.  6,  November-Deeember  1970,  p.541-549. 

26.  Shankar,  V.,  "Optimal  Estimation  and  Control  of  Linear  Time- 

Invariant  Systems  with  Uncertain  Parameters,"  Ph.D. 

Dissertation,  Stanford  University,  March  1980. 

27.  Tlroshi,  1.,  and  Elliot,  J.A.,  "Explicit  node)  Following  Control 

Scheme  Incorporating  Integral  Feedback,"  Journal  g£  Aircraft, 
Vol.  11,  No.  6,  June  1974,  p.364-366. 

28.  Trankle,  T.L.,  and  Bryson,  A.E.,  "Control  Logic  to  Track  Outputs  of 

a  Command  Cenerator,"  Journal  gi  Guidance  and  Control ,  Vol.  1, 
No.  2,  March-April  1978,  p.130-135. 

29.  Tyler.  J.S.  Jr.,  "The  Characteristics  of  Model-Following  Systems  as 

Synthesized  by  Optimal  Control,"  IEEE  Transactions  gg  Automatic 
Control .  Vol.  AC-9,  No.  4,  October  1964,  p.485-498. 

30.  Vinkler,  A.P.,  "Optimal  Controller  Design  Methods  for  Linear 

Systems  with  Uncertain  Parameters  -  Development,  Evaluation,  and 
Comparison,"  Ph.O.  Thesis,  California  Institute  of  Technology, 
April  2.  1979. 

31.  Uinsor,  C.A.,  and  Roy,  R.J.,  "The  Application  of  Specific  Optimal 

Control  to  the  Oesign  of  Desensitized  Model  Following  Control 
Systems,"  IEEE  Transactions  on  Automatic  Control .  Vol.  AC-15, 

No.  3,  June  1970,  p.326-333. 

32.  Yore,  E.E.,  "Optimal  Decoupling  Control."  Preprints  of  the  Joint 

Automatic  Control  Conference,  Ann  Arbor,  Michigan,  June  1968, 
p. 327-336. 


